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1 . 0 INTRODUCTION 



This Users Guide explains how to use the new version of the 
Naval Postgraduate School Random Number Generator Package [1], called 
LLRANDOMII. It is an outgrowth of the old LLRANDOM package [1] which 
has been the subject of several modifications and enhancements during 
the period 1978 to 1981, and which was based on the work of Lewis, 
Goodman & Miller [2]. 

This document does not go into the details of how the different 
methods operate internally, but concentrates on helping the reader to 
be able to use the package. The details are discussed in a separate 
report [3]; other references are given there. However, the philosophy 
behind the package is to provide arrays of uniform and non-uniform 
random numbers which have ; 

(i) known and documented statistical properties [6]: 

(ii) are efficiently computed; and 

(iii) are reliably computed. 

With point (iii) in mind, nothing in LLRANDOM was changed unless it 
significantly enhanced the capabilities of the package. 

This LLRANDOMII package is mostly written in IBM/360 assembly 
language but uses two subroutines written in standard IBM FORTRAN. 

Since IBM/360 Assembly language is upward compatible, the program will 
run on IBM/370 machines and on the IBM 3033. Testing was done on an 
IBM/360/65 and on an IBM 3033. All the programs are independent modules 
which follow all IBM standard linkage conventions and therefore can be 
called from all IBM high level languages which use the standard linkage 
conventions (FORTRAN, PL/I, etc.). 

The entire package is single precision oriented, so all the 
parameters are either full word integers or short precision floating 
point . Beside integer and floating point pseudo-random deviates the 
package can generate Normal, Exponential, Gamma, Cauchy, Poisson, and 
Geometric random deviates, or the order statistics associated with 
these random deviates. 



1.1 



2.0 



DIFFERENCES FROM THE PREVIOUS PACKAGE 



The new version of LLRANDOM (LLRANDOMII) is not compatible with 
the old version. Therefore users who may want to convert to the new 
package should change their subroutine name and argument list since 
both have changed from LLRANDOM to LLRANDOMII (subroutine names in 
LLRANDOMII are all different from corresponding subroutine names in 
LLRANDOM) . 

The new package has a new parameter (except on integer uniform 
generators), the fifth one on the argument list, which must be included. 
This parameter must be either integer 0 to indicate that the array of 
random numbers is not to be sorted or integer 1 to request that the 
array of random numbers be sorted in ascending order. Table 1 contains 
all the subroutine names for all generators available with the argument 
list included. 

The 4th argument in all generators indicates the multiplier to be 
used in the congruential generator. The table of multipliers currently 
has two choices (1 = 16807, 2 = 397204094). There are two versions of 
each generator, the second providing deviates based on a self-shuf f led 
version of the basic random number generator [3]. 

Using the long multiplier results in "better" random, numbers at 
the expense of an increase in execution time (up to twice as much in 
an IBM/360. No difference in time found in the IBM 3033. The shuffling 
feature also improves the quality of the "uniform" random numbers gen- 
erated, and hopefully improves the quality of the resulting non-uniform 
numbers generated when that is the case. Shuffling adds only a very 
small amount of time when used and allows the generation of sequences 
of uniforms of very long cycle. Given the availability of the long 
m.ultiplier and the shuffling feature it may have been wise to drop the 
short multiplier. However this multiplier, originally proposed by 
Lewis, Goodman & Miller [2], is widely used and is only marginally de- 
fective. Thus it can be used to rerun old simulations or to obtain in 
this package two different pseudo-random number streams. 

The new package does not require the initial call to the OVFLOt*’ 
subroutine to establish the interrupt environment. This is handled 
internally now by the LLRANDOMII programs, and represents a significant 
im.provement . 

New generators have been added for Geometric and Poisson deviates, 
enhancing the Normal, Gamma, Exponential and Cauchy capabilities of 
LLRANDOM and its extensions [5]. A new faster algorithm is used to 
generate Gamma deviates when a > 1. 

2.1 



SUBROUTINE NAMES IN LLRANDOMII PACKAGE 
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3.0 RELATIONSHIP TO IMSL (EDITION 8) , CHAPTER G 

The routines in the International Mathematics and Statistics 
Libraries (IMSL) Edition 8 (designed by P.A.W. Lewis & J. Gentle) 
implement substantially the same algorithms as are implemented in 
LLRANDOMII. However the implementation in the IMSL Library is meant 
first of all to be portable across the computing systems which IMSL 
supports, and secondly is oriented toward generation of random numbers 
one at a time. Thus for efficiency the LLRANDOMII package should be 
used. 

Timings for the LLRANDOMII routines are given in Table 2 below 
for the IBM 3033. Some timings are given for IMSL subroutines. In 
all cases timings per deviate are based on generation of arrays of 

10.000 deviates. Timings for LLRANDOM were given in [1]. LLRANDOMII 
on an IBM 3033 is roughly eight times as fast as LLRANDOMII on the 
IBM 360/67. 

TABLE 2 

TIMING COMPARISONS AND SUBROUTINE EQUIVALENCES 



LRANDOMII 



IMSL 



VARIATE 


IBM 3033 




IBM 3033 




INTEGER UNIFORM 


LINT 


(2) 




GGUD 


(16) 




REAL UNIFORM 


LRND 


(3) 




GGUBS 


(9) 




NORMAL 


LNORti 


(7) 




GGNML 


(44) 




EXPONENTIAL 


LEXPN 


(7) 




GGEXN 


(23) 








P 

II 

• 

00 


(54) 


1 


ra= . 8 


(91) 


GAMJIA (CHI-SQUARE) 


LGAMA 


U=2.5 


(18) 


GGAMR j 


^a=2.5 


(105) 






“N 

II 

00 


(5) 




II 

• 

00 


(21) 


GEOMETRIC 


LGEOM 


"^p=.98 


(16) 


GGEOT j 


^p=.98 


(21) 






II 

• 

00 


(5) 




II 

00 


(29) 


POISSON 


LPOIS 


^P=20. 


(7) 


GGPON ^ 


'■p=20. 


(323) 


CAUCHY 


LCCHY 




(7) 


GGCAY 




(34) 



— 6 

Time in microseconds (10 secs) per deviate based on generation 
of arrays of 10,000 deviates at a time is given in parenthesis. 
The Assembly Language LLRANDOMII routines are usually 3 to eiaht 
times as fast as the FORTRAN IMSL subroutines. 

For times to generate deviates one at a time add approximately 
7 microseconds to the values for LRANDOMII subroutines, 12 to 20 
microseconds to IMSL, subroutines. 

Use of the shuffled generators adds about 5% to the times for 
LRANDOMII subroutines. 
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4.0 USING THE PACKAGE 



For most of the uses initialization of the "seed" IX is the 
only thing required in addition to the call to the appropriate sub- 
routine. There are some cases though when the user may want to 
generate several streams of numbers and then he must be aware of 
some rules to follow if he is concerned about being able to repeat 
the exact sequence of numbers later on. The "seed" is the starting 
value for the (integer) uniform random number generator and it must 
be initialized to any value between 1 and 2^^-l. It is updated 
automatically on every call to the generators so that next time the 
subroutine is called the sequence of numbers produced will continue 
at the point it ended last time . Normally, the user should never 
change the seed in the course of the program. If the user wants dif- 
ferent streams of numbers at different points in his program, then he 
should use different seed names and initialize them with different values 
for the different calling points. 

In the new version of LLRANDOM all the generators are independent 
of each other since each generator is now an independent module with 
its own separate storage. However, using the same generator at different 
points in the same program will still cause some interaction in the 
case of "shuffled" generators. This is because only one copy of the 
generator is loaded by the link-editor and the shuffling table there- 
fore becomes the same for the different calling points. Then, in this 
case, if repeatability is important the user must make sure that the 
"seeds" and the number of deviates requested at every calling point 
remain the same. Section 4.1 shows an example of a case where the 
sequences do not match in different runs of the program due to a "small" 
modification made to it. As indicated before this happens only when 
using shuffled generators. 

4.1 EXAMPLES 

a) To generate 10,000 geometric deviates with parameter p=0.80, 
using the short multiplier, no shuffling, no sorting. The 
Geometric generator produces deviates X according to the 
following form of the Geometric distribution; 



f (k)=P{X=k}=p^(l-p) , 



k=0,l 



9***9 



which has moments 



E(X)=p/(l-p) , 



var (x)=p/(l-p) ^ . 



The FORTRAN program may look like this 
REAL *4 A(IOOOO) , P 
Integer IX, N 
IX=15987 
N=10000 

P=.80 

• 

Call LGEOM (IX , A,N, 1 , 0 , P) 

# 

End. 

Thus there is no sorting and the Multiplier 1 is used. 

b) Generate two batches of ten samples each of Poisson deviates. 
Batch No. 1 with samples of size Nl=5000 each and Batch No. 2 
with samples of size N2=3500 numbers. Both batches are 
from the Poisson distribution with X = 3.7, using multiplier 
2, shuffled and no sorting. 

The Poisson distribution has the form 



f (k) =P{x=k}=X^e“^/ (k! ) 



with moments 



E(x)=X, 



var (x) =X . 



4 . 2 



c 



c 



100 



The FORTRAN program may be set up like this: 
Real *4 Al (5000) ,A2 (4000), AL 
Integer IXl, 1X2, Nl, N2 
1X1=14872 
1X2=932789 
Nl=5000 
N2=3500 
AL=3 . 7 



Do 100 ISAMPL=1,10 



Calling sequence for samples of Batch No. 1 
CALL SPOIS (IX1,A1,N1,2,0,AL) 

Calling sequence for samples of Batch No. 2 
CALL SPOIS (IX2,A2,N2,2,0,AL) 

CONTINUE 

END 

This example shows one case where the same shuffled 
generator is used in different points in the sane program, 
and interaction between the two calls is created. If the 
same sequences of numbers have to be repeated in a future 
run of the program, the user must make sure that 1X1, 1X2, 

Nl and N2 remain the sam.e. For example, if the programt is 
run again with N2=4000 the shuffling table (which is shared 
by the two calls) will be updated more often now in the 
second call causing the numbers from the first call to be 
different. 

As indicated before this happens only when the same sub- 
routine name is used in both calls and they are of the 
shuffled type. 



4.3 



As far as the no-shuffle type generators being able to 
regenerate the same sequences of numbers, all that is 
required is to use the same value for the "seed" in all 
different runs of the program. Different names for the 
seed variable are still necessary to avoid interaction 
between the different calls, in the case of multiple 
calls in the same user program. 

4.2 NOTE FOR PL/1 USERS 

When using the LLRANDOMII package from PL/1 programs the ap- 
propriate subroutines being used must be declared as external entries 
with data types Binary Fixed (31) and Binary Float (21) for the 
parameters to be passed that correspond, respectively, to INTEGER 
and REAL *4 in Fortran. For example, if the LGAMA generator is used 
in a PL/1 program, the following declaration must be included in the 
program; 

DCL LGAMA EXTERNAL ENTRY (BIN FIXED (31) , (*) BIN FLOAT ( 21) , 

BIN FIXED(31), BIN FIXED(31), BIN FIXED(31), BIN FL0AT(21)) 

OPTIONS (ASSEMBLER INTER) ; 

5.0 AVAILABILITY 

The LLRANDOMII package is distributed in 9 track 1600BPI m.agnetic 
tape containing two partitioned data sets. The first contains the 
source IBM Assembly and Fortran code and the second contains object 
code that can be used to construct a TXTLIB under VM/CMS or as a 
library in the link-edit step when used in batch mode. 

To obtain a copy send a scratch tape to 

Naval Postgraduate School 
Monterey, CA 93940 
Operations Research Dept. 

Attn; Professor P.A.W. Lev/is 
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